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Abstract 


This  report  describes  a  new  filtering  approach  suitable  for  demodulating  narrowband 
signals  captured  in  digital  wideband  data.  The  new  approach  is  based  on  generating  the 
filtered  and  downconverted  baseband  signals  using  time/frequency  spectral  data,  rather 
than  the  input  wideband  data.  For  communications  electronic  support  measures 
systems,  this  spectral  data  will  normally  be  generated  to  facilitate  signal  detection  and 
parameter  estimation  (i.e.  frequency,  bandwidth,  bearing,  etc.).  Hence  the  new 
approach  is  able  to  take  good  advantage  of  this  to  reduce  the  overall  number  of 
computations.  Through  comparisons  with  more  traditional  FIR  filtering  techniques,  it  is 
shown  that  the  new  spectral  filtering  approach  is  computationally  more  efficient  when 
multiple  signals  are  to  be  simultaneously  demodulated. 


Resume 


Le  present  rapport  decrit  une  nouvelle  methode  de  filtrage  qui  convient  a  la 
demodulation  de  signaux  a  bande  etroite  saisis  dans  des  donnees  numeriques  a  large 
bande.  La  nouvelle  approche  repose  sur  la  production  des  signaux  en  bande  de  base 
fibres  et  abaisses  au  moyen  de  donnees  spectrales  temps-frequence,  plutot  qu’au  moyen 
de  donnees  d’ entree  a  large  bande.  Dans  le  cas  des  systemes  de  mesures  de  soutien 
electronique  (MSE)  pour  les  communications,  ces  donnees  spectrales  sont 
normalement  generees  pour  faciliter  la  detection  des  signaux  et  T estimation  des 
parametres  (conime  la  frequence,  la  largeur  de  bande  et  le  relevement).  La  nouvelle 
approche  peut  done  en  tirer  profit  pour  reduire  le  nombre  global  de  calculs.  Par  des 
comparaisons  aux  techniques  de  filtrage  FIR  traditionnelles,  on  montre  que  la  nouvelle 
methode  de  filtrage  spectral  est  plus  efficace  sur  le  plan  des  calculs  lorsqu'il  faut 
demoduler  plusieurs  signaux  simultanement. 
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Executive  summary 


Spectrum  monitoring  and  communications  electronic  support  measures  systems  (ESM) 
require  considerable  flexibility  and  sophistication  to  cope  with  increasingly  complex 
signal  environments.  One  system  designed  to  meet  these  needs  is  the  Military  Digital 
Analysis  System  (MiDAS).  It  is  a  wideband  intercept/DF  system  which  uses  a 
combination  of  personal  computer  based  digital  signal  processing  algorithms  and  a 
modular  system  architecture  based  on  commercial-off-the-shelf  (COTS)  hardware. 
Testing  carried  out  to  date  has  shown  the  system  has  the  potential  to  meet  current  and 
future  requirements,  and  current  research  and  development  is  being  carried  out  to 
achieve  its  full  functionality. 

One  capability  of  interest  is  the  simultaneous  demodulation  of  narrowband  signals  of 
interest.  From  the  hardware  perspective,  MiDAS  is  well  suited  for  this  application  since 
a  wideband  receiver  can  be  set  up  to  search  the  frequency  band  where  targets  of  interest 
arc  expected.  The  difficulty  is  that  the  wider  the  instantaneous  frequency  coverage,  the 
greater  the  data  rate  at  the  output  of  the  receiver,  and  the  greater  the  likelihood  of 
intercepting  extraneous  interference  and  noise.  Sorting  through  all  this  data  to 
determine  the  signal  parameters  of  the  desired  signals  is  not  a  trivial  task,  especially  in 
a  rapidly  changing  signal  environment. 

A  common  approach  to  dealing  with  wideband  signal  data  is  to  divide  the  data  into 
blocks  of  N  samples  and  then  convert  each  block  of  data  to  the  frequency  domain  using 
an  FFT  processor.  The  result  is  a  series  of  spectral  snapshots  which  change  over  time  as 
the  signal  environment  changes,  i.e.  a  time/frequency  spectrum  is  produced.  This 
works  well  for  the  detection  and  estimation  of  narrowband  fixed  frequency  signals 
since  they  occupy  small  portions  of  the  spectrum  at  any  instance  in  time  and  can  be 
separated  based  on  frequency,  while  noise  is  spread  out  across  the  entire  spectrum 
weakening  its  disruptive  influence  accordingly. 

Once  signals  have  been  detected,  estimation  of  the  signal  parameters  follows  in  order  to 
determine  if  the  signal  is  of  interest  and  to  acquire  the  parameters  needed  for 
demodulation  purposes.  Before  demodulating  the  signal,  the  data  is  filtered  to  remove 
noise  and  any  other  signals  which  might  interfere  with  the  demodulation.  One  way  to 
do  this  is  to  pass  the  wideband  data  through  a  bandpass  filter  which  has  been  matched 
to  the  frequency  and  bandwidth  of  the  targeted  signal. 

Although  this  approach  works  well,  it  may  not  be  the  most  computationally  efficient 
way  given  that  filtered  data  is  already  available  in  the  form  of  the  time/frequency 
spectrum.  In  this  report,  an  approach  (called  the  spectral  filter  here)  is  proposed  which 
uses  the  time/frequency  spectral  data  to  generate  the  required  filtered  data.  The  main 
requirements  to  implement  this  approach  are:  a  50%  overlap  when  generating  the 
time/frequency  spectrum  (i.e.  the  start  of  each  data  block  used  for  a  spectral  snapshot 
overlaps  the  previous  block  by  50%);  and  the  sum  of  two  FFT  windows  overlapped  by 
50%  is  unity.  The  downside  of  these  requirements  is  that  it  doubles  the  amount  of 
processing  versus  the  case  when  no  overlap  is  used  (which  is  how  MiDAS  is  currently 
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operated),  and  it  restricts  the  choices  of  FFT  windows  that  can  be  used.  However,  if 
optimum  signal  detection  is  a  goal,  overlapping  will  be  required  as  a  matter  of  course 
[1],  and  the  restriction  on  the  choice  of  FFT  windows  is  not  a  debilitating  condition. 

The  advantage  of  the  spectral  filter  is  that  the  bulk  of  the  processing  involves  the 
generation  of  the  time/frequency  spectrum  which  can  be  done  efficiently  using  a  fast 
Fourier  transform  (FFT).  Subsequent  processing  to  extract  and  filter  the  desired  signals 
requires  very  little  additional  computation,  hence  there  is  relatively  little  difference 
whether  one  or  a  number  of  signals  are  extracted.  By  way  of  comparison,  for  normal 
time  domain  based  filters  the  processing  increases  as  a  linear  function  of  the  number  of 
signals  and  may  become  prohibitive  if  more  than  a  few  signals  are  to  be  demodulated. 

For  comparative  purposes,  two  other  filtering  approaches  were  also  introduced  for 
study,  namely,  a  single  stage  FIR  filter  and  a  multistage  FIR  filter.  MiDAS  currently 
uses  a  multistage  FIR  filter.  The  single  stage  filter  has  the  greatest  flexibility  in  terms  of 
designing  the  filter  frequency  response.  The  multistage  filter  has  considerably  less 
flexibility,  but  is  computationally  more  efficient  than  the  single  stage  filter. 

In  terms  of  designing  the  filter  frequency  response,  the  spectral  filter  falls  in  between 
the  single  stage  and  multistage  filters.  It  is  more  flexible  and  has  a  better  response  than 
the  multistage  filter,  but  is  still  relatively  inflexible  when  compared  to  the  single  stage 
filter. 

In  terms  of  computational  speed,  both  the  single  stage  and  multistage  filters  are  faster 
when  only  one  signal  is  processed.  However  when  two  signals  are  processed 
simultaneously,  the  spectral  filter  is  faster  than  the  single  stage  filter,  and  when  four 
signals  are  processed  simultaneously,  the  spectral  filter  is  also  faster  than  the  multistage 
filter. 

The  conclusion  is  that  for  typical  applications  involving  modern  wideband  ESM 
systems,  such  as  MiDAS,  where  demodulation  of  four  or  more  narrowband  signals  is 
desirable  and  computations  are  at  a  premium,  then  the  spectral  filtering  approach 
proposed  in  this  report  is  the  preferred  choice. 
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Sommaire 


Les  systemes  de  mesures  de  soutien  electronique  (MSE)  pour  les  communications  et  la 
surveillance  du  spectre  ont  besoin  d’un  perfectionnement  pousse  et  d’une  grande 
souplesse  pour  fonctionner  dans  des  conditions  de  transmission  des  signaux  de  plus  en 
plus  complexes.  Un  systeme  concu  pour  repondre  a  ces  besoins  est  le  systeme 
d’ analyse  numerique  militaire  (MiDAS).  II  s’agit  d’un  systeme  de 
radiogoniome  trie/d’ interception  a  large  bande  qui  a  recours  a  la  fois  a  des  algorithmes 
de  traitement  numerique  sur  ordinateur  personnel  et  a  une  architecture  de  systeme 
modulaire  sur  materiel  commercial  courant.  Les  essais  menes jusqu’a  maintenant 
montrent  que  le  systeme  offre  la  possibility  de  repondre  aux  besoins  actuels  et  a  venir, 
et  des  travaux  de  recherche  et  developpement  sont  menes  en  vue  de  la  realisation 
integrate  de  cette  fonctionnalite. 

Une  capacite  souhaitable  est  la  demodulation  simultanee  de  signaux  a  bande  etroite 
presentant  un  interet.  Du  point  de  vue  du  materiel,  le  systeme  MiDAS  convient  bien  a 
cette  application,  du  fait  qu’un  recepteur  large  bande  peut  etre  configure  pour  chercher 
la  bande  de  frequences  dans  laquelle  on  s’ attend  a  des  cibles  presentant  un  interet.  La 
difficulty,  e’est  que,  plus  la  frequence  de  couverture  instantanee  est  large,  plus  le  debit 
est  eleve  a  la  sortie  du  recepteur  et  plus  il  y  a  de  chances  d'intercepter  du  bruit  et  du 
brouillage  etrangers.  Le  triage  des  donnees  en  vue  de  determiner  les  parametres  des 
signaux  utiles  est  une  tache  non  negligeable,  surtout  lorsque  les  signaux  changent 
rapidement. 

Une  approche  commune  au  traitement  des  donnees  relatives  aux  signaux  a  large  bande 
consiste  a  diviser  les  donnees  en  blocs  de  N  echantillons,  puis  a  convertir  chaque  bloc 
de  donnees  dans  le  domaine  frequence  au  moyen  d’un  processeur  LET.  II  en  resulte  une 
serie  d’instantanes  spectraux,  qui  changent  au  fur  et  a  mesure  de  revolution  des 
conditions  des  signaux,  e’est-a-dire  qu’un  spectre  temps-frequence  est  produit.  Cela 
convient  bien  a  la  detection  et  a  l’estimation  des  signaux  a  frequence  fixe  a  bande 
etroite,  du  fait  qu’ils  occupent  de  faibles  portions  du  spectre  a  n'importe  quel  moment 
et  qu’ils  peuvent  etre  separes  en  fonction  de  la  frequence,  tandis  que  le  bruit  est  etale 
dans  1’ ensemble  du  spectre,  ce  qui  reduit  son  influence  perturbatrice. 

Une  fois  les  signaux  detectes,  on  procede  a  une  estimation  de  leurs  parametres  pour 
determiner  s’ils  presentent  un  interet  et  pour  acquerir  les  parametres  necessaries  aux 
fins  de  demodulation.  Avant  que  les  signaux  ne  soient  demodules,  on  filtre  les  donnees 
pour  en  retirer  le  bruit  et  tout  autre  signal  susceptible  de  perturber  la  demodulation.  A 
cette  fin,  on  peut  faire  passer  les  donnees  a  large  bande  dans  un  filtre  passe-bande  regie 
a  la  frequence  et  a  la  largeur  de  bande  du  signal  cible. 

Meme  si  cette  approche  donne  de  bons  resultats,  elle  n’est  pas  necessairement  la  plus 
efficace  sur  le  plan  des  calculs,  puisque  les  donnees  filtrees  sont  deja  disponibles  sous 
la  forme  de  spectre  temps-frequence.  Dans  le  present  rapport,  on  propose  une  approche 
(appelee  filtrage  spectral)  qui  fait  appel  aux  donnees  spectrales  temps-frequence  pour 
general-  les  donnees  filtrees  requises.  Les  principals  exigences  de  la  mise  en  oeuvre  de 


DRDC  Ottawa  TM  2005-198 
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cette  approche  sont  les  suivantes  :  chevauchement  de  50%  au  moment  de  generer  le 
spectre  temps-frequence  (c’est-a-dire  que  le  debut  de  chaque  bloc  de  donnees  utilise 
pour  l’obtention  d’un  instantane  spectral  chevauche  de  50%  le  bloc  precedent);  et  la 
somme  de  deux  fenetres  FFT  chevauchees  de  50%  est  l’unite.  L’ inconvenient  de  ces 
exigences,  c’est  qu’elles  doublent  le  volume  de  traitement  par  rapport  au  cas  oil  il  n’y  a 
aucun  chevauchement(c’est  de  cette  facon  que  le  systeme  MiDAS  fonctionne  a  l’heure 
actuelle),  et  elles  limitent  la  selection  de  fenetres  FFT  dont  on  peut  se  servir.  Toutefois, 
si  Fobjectif  est  la  detection  optimale  de  signaux,  le  chevauchement  constitue  une 
exigence  tout  a  fait  normale  [citer  le  document  de  reference  1],  et  la  restriction  imposee 
a  l’egard  de  la  selection  des  fenetres  FFT  ne  constitue  pas  une  condition  debilitante. 

L’avantage  du  filtre  spectral,  c’est  que  la  majeure  partic  du  traitement  comporte  la 
generation  du  spectre  temps-frequence,  ce  qui  peut  se  faire  efficacement  a  l’aide  d’une 
transformee  de  Fourier  rapide  (FFT).  Le  traitement  subsequent  en  vue  de  Textraction  et 
du  filtrage  des  signaux  utiles  requiert  tres  peu  d’operations  additionnelles,  ce  qui  veut 
dire  qu’il  n’y  a  pas  tellement  de  difference  a  extraire  un  certain  nombre  de  signaux  ou 
un  seul.  Par  comparaison,  dans  le  cas  des  fibres  a  reponse  temporelle  normaux,  le 
traitement  augmente  en  fonction  lineaire  du  nombre  de  signaux  et  risque  de  devenir 
prohibitif  s’il  faut  demoduler  plus  que  quelques  signaux. 

Aux  fins  de  comparaison,  deux  autres  approches  de  filtrage  ont  egalement  ete  etudiees, 
a  savoir  un  filtre  FIR  a  un  seul  etage  et  un  filtre  FIR  a  plusieurs  etages.  Le  systeme 
MiDAS  utilise  actuellement  un  filtre  FIR  a  plusieurs  etages.  C’est  le  filtre  a  un  seul 
etage  qui  offre  la  plus  grande  souplesse  en  ce  qui  concerne  la  conception  de  la  reponse 
en  frequence  du  filtre.  Le  filtre  a  plusieurs  etages  est  beaucoup  moins  souple,  mais  il  est 
plus  efficace  sur  le  plan  des  calculs  que  le  filtre  a  un  seul  etage. 

Pour  ce  qui  est  de  la  reponse  en  frequence  souhaitee,  le  filtre  spectral  se  situe  entre  le 
filtre  a  un  seul  etage  et  le  filtre  a  plusieurs  etages.  Il  est  plus  souple  et  a  une  meilleure 
reponse  que  le  filtre  a  plusieurs  etages,  mais  il  reste  relativement  peu  souple  par  rapport 
au  filtre  a  un  seul  etage. 

En  ce  qui  concerne  la  vitesse  de  calcul,  le  filtre  a  un  seul  etage  et  le  filtre  a  plusieurs 
etages  sont  plus  rapides  lorsqu’un  seul  signal  est  traite.  Par  contre,  lorsque  deux 
signaux  sont  traites  simultanement,  le  filtre  spectral  est  plus  rapide  que  le  filtre  a  un 
seul  etage  et,  lorsque  quatre  signaux  sont  traites  simultanement,  le  filtre  spectral  est 
aussi  plus  rapide  que  le  filtre  a  plusieurs  etages. 

La  conclusion  est  que,  pour  les  applications  typiques  comportant  des  systemes  MSE  a 
large  bande  modernes,  comme  le  systeme  MiDAS,  ou  la  demodulation  de  quatre 
signaux  a  bande  etroite  ou  plus  est  souhaitable  et  ou  la  capacite  de  calcul  est  limitee,  le 
filtrage  spectral  propose  dans  le  present  rapport  represente  T  approche  de  choix. 


Read,  W.J.L.  2005.  A  Time/Frequency  Spectral  Based  Approach  to  Filtering.  DRDC  Ottawa 
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VI 


DRDC  Ottawa  TM  2005-198 


Table  of  contents 


Abstract .  i 

Resume .  i 

Executive  summary .  iii 

Sommaire .  v 

Table  of  contents  .  vii 

List  of  figures . viii 

List  of  tables .  ix 

1.  Introduction .  1 

2.  Liltering  using  Time/Lrequency  Spectral  Data .  3 

2.1  Time/Lrequency  Spectral  Estimation .  3 

2.2  Generating  the  Time  Domain  Liltered  Data .  4 

2.3  Processing  Requirements .  10 

3.  Temporal  Liltering .  14 

3.1  Single  Stage  Lilter .  14 

3.2  Multistage  Liltering .  15 

3.2.1  Halfband  Lilter .  15 

3.2.2  Composite  Lilter .  19 

4.  Comparative  Results .  21 

4.1  Lrequency  Response .  21 

4.2  Computational  Comparisons .  24 

5.  Conclusions .  27 

References .  29 

Annex .  31 

DRDC  Ottawa  TM  2005-198  vii 


List  of  figures 


1  Time/frequency  spectrum  example .  4 

2  Impulse  noise  filtering  example .  7 

3  Filter  response  as  a  function  of  filter  length .  21 

4  Filter  gain  at  a  normalized  frequency  of  ±0.75  as  a  function  of  filter  length  ....  22 

5  Frequency  response  of  three  different  filtering  approaches .  23 

6  Computation  ratio  as  a  function  of  filter  length .  25 

7  Computation  ratio  as  a  function  of  the  oversampling  factor .  25 

8  Parameter  curves  (/3  versus  n) .  27 


viii 


DRDC  Ottawa  TM  2005-198 


List  of  tables 


1  Adjacent  channel  filter  response .  9 

DRDC  Ottawa  TM  2005-198  ix 


This  page  intentionally  left  blank. 


x 


DRDC  Ottawa  TM  2005-198 


1.  Introduction 


Spectrum  monitoring  and  communications  electronic  support  measures  systems  require 
considerable  flexibility  and  sophistication  to  cope  with  increasingly  complex  signal 
environments.  One  system  designed  to  meet  these  needs  is  the  Military  Digital  Analysis 
System  (MiDAS).  This  is  a  wideband  intercept/DF  system  which  uses  a  combination  of 
personal  computer  based  digital  signal  processing  algorithms  and  a  modular  system 
architecture  based  on  commercial-off-the-shelf  (COTS)  hardware.  Testing  carried  out 
to  date  has  shown  the  system  has  the  potential  to  meet  current  and  future  requirements, 
and  current  research  and  development  is  being  carried  out  to  extend  its  capabilities. 

One  capability  of  interest  is  the  simultaneous  demodulation  of  narrowband  signals  of 
interest.  From  the  hardware  perspective,  MiDAS  is  well  suited  for  this  application  since 
a  wideband  receiver  can  be  set  up  to  cover  large  portions  of  the  frequency  band  where 
targets  of  interest  arc  expected  to  occur.  The  difficulty  is  that  the  wider  the  frequency 
coverage,  the  greater  the  data  rate  at  the  output  of  the  receiver,  and  the  greater  the 
likelihood  of  intercepting  extraneous  interference  and  noise.  Sorting  through  all  this 
data  to  determine  the  parameters  of  the  desired  signals  is  not  a  trivial  task,  especially  in 
a  rapidly  changing  signal  environment. 

A  common  approach  to  dealing  with  wideband  data  is  to  divide  the  data  into  blocks  of 
N  samples  and  then  transform  each  block  of  data  samples  to  the  frequency  domain 
using  a  Fast  Fourier  transform  (FFT).  The  result  is  a  series  of  spectral  snapshots  which 
change  over  time  as  the  signal  environment  changes,  i.e.  a  time/frequency 
representation  of  the  spectrum  is  produced.  This  works  well  for  the  detection  and 
estimation  of  narrowband  fixed  frequency  signals  since  they  occupy  small  portions  of 
the  spectrum  at  any  instance  in  time  and  can  be  separated  based  on  frequency,  while 
noise  is  spread  out  across  the  entire  spectrum  weakening  its  disruptive  influence 
accordingly. 

Once  signals  have  been  detected,  estimation  of  the  signal  parameters  follows  in  order  to 
determine  if  the  signal  is  of  interest  and  to  acquire  the  parameters  needed  for 
demodulation.  Before  demodulating  the  signal,  the  data  is  filtered  to  remove  noise  and 
any  other  signals  which  might  interfere  with  the  demodulation.  One  way  to  do  this  is  to 
pass  the  wideband  data  through  a  bandpass  filter  which  has  been  matched  to  the 
frequency  and  bandwidth  of  the  targeted  signal. 

Although  this  approach  works  well,  it  may  not  be  the  most  computationally  efficient 
way  to  do  it  given  that  filtered  data  is  already  available  in  the  form  of  the  time/frequency 
spectrum.  In  this  report,  an  approach  is  proposed  which  uses  the  time/frequency  data  to 
generate  the  required  filtered  data.  The  main  requirements  to  implement  this  approach 
arc:  a  50%  overlap  of  the  data  blocks  when  generating  the  time/frequency  spectrum 
(i.e.  the  start  of  each  data  block  used  for  a  spectral  snapshot  overlaps  the  previous  block 
by  50%);  and  the  sum  of  two  FFT  windows  overlapped  by  50%  is  unity. 

In  the  rest  of  this  report.  Section  2  introduces  a  new  approach  to  filtering  using  the 
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time/frequncy  data  and  discusses  the  processing  requirements.  Section  3  introduces 
single  stage  FIR  filtering  and  multistage  FIR  filtering  approaches  for  the  purposes  of 
comparison  with  the  new  approach.  This  section  also  discusses  the  processing 
requirements  of  the  FIR  approaches.  Section  4  provides  comparisons  between  the  new 
filtering  approach  and  the  FIR  approaches,  both  in  terms  of  filter  frequency  response 
and  processing  requirements.  Finally,  Section  5  provides  the  concluding  remarks. 
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2.  Filtering  using  Time/Frequency  Spectral  Data 

2.1  Time/Frequency  Spectral  Estimation 


A  block  of  wideband  data  is  represented  here  as 

X  =  [xi,X2,  ...,XN]T 


(1) 


where  the  superscript  T  represents  transpose  and  x  \ ,  x<> , . . . ,  xjv  represent  the  complex 
sequential  outputs  from  the  wideband  receiver  (some  pre-processing,  such  as  frequency 
shifting,  anti-aliasing  filtering,  and  conversion  from  real  to  complex  formatted  data, 
may  have  already  taken  place).  To  produce  a  time/frequency  spectrum,  the  data  is 
rearranged  into  columns  defined  by  the  M  x  K  matrix 


W \X\ 

WlMXk+l 

WlX2k+ 1 

■  ■  WiXN-M+1 

W2X2 

W2Xk+2 

W2X2k+2 

■  ■  W2XN-M+2 

WMXM 

WMXk+M 

WMX2k+M  ' 

WMXN 

where  M  is  the  number  of  frequency  bins  or  channels  desired  for  each  spectral 
snapshot  and  is  assumed  to  be  a  power  of  two  (for  computational  efficiency  of  the 
FFT),  the  weighting  coefficients  w\ ,u>2,  ■■■,  wm  represent  the  FFT  window,  M  —  k 
represents  the  overlap  (k  =  M  for  no  overlap),  and  the  number  of  columns  is  given  by 


K 


N  -  M 
k 


+  1 


(3) 


For  the  approach  proposed  in  this  report,  a  50%  overlap  is  used  which  corresponds  to 
k  =  M/2  where  M  is  even- valued.  The  number  of  columns,  then  becomes 


K 


2N 

~M 


(4) 


For  convenience,  it  is  assumed  that  N  is  chosen  to  exactly  satisfy  the  expression 


K 


2N 

~M 


(5) 


where  K  is  a  positive  even- valued  integer.  Although  other  values  of  N  arc  permitted, 
the  extra  data,  xM(k+1y2+i,  ...,  xjv  and  the  last  column  of  Y  (if  odd-numbered),  would 
effectively  be  ignored. 


Once  Y  has  been  created,  the  FFT  is  applied  to  each  column  to  convert  the  column 
time  domain  data  to  frequency  domain  data.  The  resultant  output  can  be  represented  as 


•Sn 

s  12  • 

■  ■  SiK 

S21 

■S22 

S2K 

SMI 

S  M2  ■ 

■  ■  SMK 
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An  example  of  the  time/frequency  spectrum  generated  using  the  approach  described 
here  is  shown  in  Figure  1.  Signal  data  was  collected  using  MiDAS  over  the  27.5  to  32.5 
MHz  band  in  a  laboratory  environment.  Three  narrowband  signals  are  featured  with  the 
third  switching  on  at  t  =  17  mS.  The  power  measurements  arc  made  relative  to  the 
median  noise  power  level  as  measured  across  the  frequency  band. 
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Figure  1:  An  example  of  a  time/frequency  spectrum  produced  from  data  featuring  two  narrowband  signals  plus  a 

third  narrowband  signal  switching  on  at  t  =  17  mS. 


2.2  Generating  the  Time  Domain  Filtered  Data 


Given  that  the  signal  of  interest  has  a  bandwidth  that  is  small  compared  to  the  receiver 
bandwidth,  then  filtering  and  decimation  of  the  data  can  be  implemented  by  limiting 
further  processing  of  the  data  to  the  rows  of  S  which  contain  the  targeted  signal 
information.  The  other  rows,  which  contain  noise  and  interference  (including  other 
signals),  are  excluded.  Hence,  denoting  the  index  numbers  of  the  signal  rows  as 
i i,  +  ii+2, ...,  12,  then  a  filtered  spectral  signal  matrix  can  be  formed  which  is  a 

subset  of  S  and  is  given  by 


Sm 


■Siil  Si\2  *  *  *  &i\K 

Si2 1  Si2  2  •  •  •  Si2K 


(V) 


where  1  <  i\  <  <  M.  An  additional  restriction  on  ii  and  i-i  is  that  the  equation 

D  =  M/(z2  —  i\  +  1)  yields  an  integer  value  for  D  where  D  is  the  decimation  rate  (the 
ratio  of  the  input  sampling  rate  to  the  output  sampling  rate  of  the  filter). 


For  demodulation  purposes,  the  time  domain  data  is  required  so  that  Sm  must  be 
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converted  accordingly.  This  is  done  in  two  steps.  The  first  step  is  to  perform  an  inverse 
FFT  (without  using  any  FFT  window)  on  the  columns  of  Sm  yielding  the  matrix 


Z 


211  212  ' ' '  2i K 


I  Z  M  I  Z  M  ry  •  •  •  g  AT  TV" 

L  D1  D2  ~DK 


(B) 


where  jj  =  i-2  —  i\  +  1.  The  columns  of  Z  represent  filtered  and  decimated  versions  of 
the  columns  of  Y.  In  the  ideal  case,  where  the  original  data  is  noise  and  interference 
free  and  the  spectral  content  of  the  signal-of-interest  is  strictly  limited  to  rows  i]  to  <9 
of  the  matrix  S,  then  the  elements  of  Z  can  be  related  back  to  the  elements  of  Y.  For 
example,  the  first  column  of  Y  is  the  inverse  Fourier  transform  of  the  first  column  of  S 
and  is  given  by 


wnxn 


1 

M 


M 

sfelej27r(fe-1)(n-1)/M 

k= 1 


(9) 


for  n  =  1, ...,  M.  Similarly,  the  first  column  of  Z  is  the  inverse  Fourier  transform  of  the 
first  column  of  Sm  which  is  given  by 


2m  1  — 


D 

M 


M 

D 

E 

i— 1 


s(i+h-l)le 


(10) 


for  m  =  1. ....  jj.  Since  the  signal  content  is  limited  to  rows  i\  to  %2  of  the  matrix  S, 
then  Ski  =  0  for  k  <  i\  or  k  >  i 2.  Given  this,  and  replacing  the  index  i  by  k  —  ii  +  1, 
then 


2ml  — 


D 

M 


M 

E 

k= 1 


®fel® 


j'27r(fc— ii)D(m—  1)  /M 


M 


=  Pe-j2,r^1“1^(m-i)/M 

M  ^ 


Sfeie 


j2n(k—l)D{m—l)/M 


(ID 


k=\ 


Comparing  the  expression  inside  the  large  parentheses  to  (9)  yields 


2m  1  =  Dpm  1W(D  m—  D+l)x(Dm-D+l)  (12) 

where  p  =  e-i27r(*i-t)£>/M  Extending  this  idea  to  all  the  columns  of  Z  then,  in  the 
special  noise  and  interference  free  case  ,  Z  can  be  rewritten  as 


=  D 


W \X\ 

W(D+ 1)  PX{D+ 1) 
W(2D+1)  P2x{2D+l) 


miX(M+i) 

W(D+ 1)  Px(M+d+ i) 
W(2D+1)  P2x(M+2D+l ) 


•  •  •  WiXN-M+1 

W(D+1)  PX(N-M+D+1) 

■  ■  ■  W(2D+l)  P  x(N-M+2D+l) 


W(M-D+ 1)  PD  X(M-D+ 1)  W(M-D+1)  FD 


W(M-D+1)  P0  ^(Af-D+l) 


(13) 
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The  second  step  in  deriving  the  time  domain  data  is  to  extract  the  required  data  from  Z. 
Using  (13),  it  is  obvious  that  the  values  of  x\,xd+i,X2d+i,  are  easily  estimated 
from  the  elements  of  Z  since  p  and  the  weighting  coefficients,  w\  ,  wn+  \ ,  W2d+\ , ..., 
arc  all  known.  The  most  obvious  method  is  by  directly  compensating  for  the  effects  of 
the  various  multiplicative  factors  according  to 


Xk% 


Zmn 

Dwk i  Pm_1 


(14) 


where  k\  =  D(m  —  1)  +  1  and  k-2  =  \M(n  —  1)  +  D(m  —  1)  +  1. 


Unfortunately  this  simple  approach  has  problems  when,  for  example,  effects  such  as 
impulsive  noise  arc  considered.  Since  the  filtering  process  described  here  is  additive 
(the  result  of  combining  several  signals  and  noise  followed  by  filtering  is  the  same  as 
individually  filtering  the  signals  and  noise  then  combining  the  results),  the  effect  of 
noise  can  be  considered  independently  from  the  effect  of  the  desired  signal.  Figure  2a 
shows  the  amplitude  of  an  impulsive  noise  signal  where  the  200th  value  is  one  and  all 
the  other  values  arc  zero.  The  noise  was  converted  to  the  frequency  domain  using  the 
Fourier  transform  with  a  Hanning  window.  Filtering  was  then  performed  assuming  the 
signal-of-interest  occupied  bins  i  \  =  65  to  12  =  72  of  the  spectral  snapshot  as  shown  in 
Figure  2b.  Using  the  inverse  Fourier  transform  to  produce  time  domain  data  yielded  the 
result  shown  in  Figure  2c.  The  interpolated  values  shown  in  the  figure  were  determined 
by  performing  zero  padding  before  applying  the  inverse  Fourier  transform.  Finally, 
Figure  2d  shows  the  result  of  correcting  for  the  Hanning  window  through  simple 
division  (e.g.  as  done  in  (14)).  The  main  observation  is  that  the  noise  values  near  the 
start  and  end  of  the  estimated  data  sequence  shown  in  Figure  2d  arc  unrealistically  large 
(i.e.  the  output  data  values  near  the  start  and  end  exceed  10000  while  the  largest  value 
in  the  input  data  never  exceeded  1).  This  is  due  to  the  fact  that  the  corresponding 
window  values  arc  very  small,  so  division  using  these  values  greatly  magnifies  the 
corresponding  noise  values. 


The  noise  enhancement  effect  can  be  quantified  in  terms  of  the  expected  change  in  the 
overall  SNR.  Noting  that  all  the  estimated  data  values  x\,xd+\,X2d+Ii  ■  •••  are 
affected,  then  the  expected  change  in  the  SNR  will  be  given  by 


A  SNR 


—10  log 


Ry _ I _ \ 

M  A-'  wf,  n  n  , . ,  / 
k= 1  (kD-D+1)  J 


(15) 


For  most  commonly  used  windows  (e.g.,  Hamming,  Hanning,  Kaiser,  Blackman,  etc.), 
the  decrease  in  SNR  is  on  the  order  of  ten’s  of  dB.  The  problem  can  be  avoided  by  not 
using  any  windowing  (or  equivalently,  letting  Wk  =  1  for  all  values  of  k),  however,  the 
sidelobe  response  of  the  filtering  process  will  be  adversely  affected. 


A  better  approach  is  to  take  advantage  of  the  fact  that  each  value  of  a;  &  for  k  =  1, ...,  N 
appeal's  twice  in  Y  except  for  the  first  and  last  M / 2  values.  That  is,  due  to  the 
overlapping  nature  of  the  data  samples,  a  value  of  X}-  appealing  in  the  latter  half  of  a 
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1.5 


Impulse  Noise  HanninaWindow 
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Figure  2:  An  example  of  filtering  impulse  noise  showing  (a)  the  input  impulse  noise  data  and  the  Hanning  window  to 
be  applied,  (b)  the  Fourier  transform  of  the  windowed  noise  (blue  line)  and  the  eight  bins  selected  to  implement  the 
filtering  (black  line),  (c)  the  inverse  Fourier  transform  of  the  eight  selected  bins  (black  lines),  as  well  as  the 
interpolated  values  (blue  line),  (d)  and  the  effect  of  correcting  for  the  Hanning  window. 
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column  Y  will  also  appeal-  in  the  first  half  of  the  next  column.  This  overlap  carries  over 
into  Z,  so  using  (13)  to  define  an  element  of  Z  as 

A?in  —  Dwkl  p  Xfc2  (16) 

where  the  row  index  m  >  M/2,  and  the  indices  k\  and  kp  ai-e  defined  as  before,  then 
the  coiTesponding  element  of  Z  containing  the  same  value  of  xk2  is  found  in  the  next 
column  and  is  defined  as 

2(rn_M)(n+ 1)  =  p^-'-ib)  Xk2  =  (-1  xk2 

(17) 

As  an  example,  consider  the  case  where  N  >  128,  M  =  32,  and  D  =  8.  For  the  first 
four  columns  of  Z,  the  elements  would  be  defined  by 


+Dw i  x\ 

j Dwg  Xg 

—Dwn  x±7  +Dw±  x\7 

+jDw2 5  X25  -jDw 9  X25 

-Dwn  %33  +Dwi  x33 

+j  D'U)25  X41  ~jDwg  X41 

—Dwn  X4g  +Dw\  X49 

+jDw 25  X57  -jDw 9  X57 

-Dwn  X(& 

+jDw 25  X73 


(18) 


where  p  =  — j,  and  the  columns  have  been  staggered  here  in  order  to  line  up  related 
values  of  Xk- 


Inspecting  the  two  elements  defined  in  (16)  and  (17),  a  useful  quantity  is  given  by  the 
sum 


Ann  +  (  l)1  Z(m-|§)(n+ 1)  —  D(Wk1  +  W(fel_M))  fj  Xfc2 

In  the  presence  of  noise,  this  expression  becomes  an  estimate.  Rearranging,  the 
corresponding  estimate  for  xp.  is  given  by 

„  Ann  T  (  l)*1 

Xk2  D{wkl+w{ki_M))pm-1 


(19) 


(20) 


which  is  similar  in  form  to  (14),  but  the  problem  term  wkl  is  replaced  by 
wkl  +  W(ki_My  For  most  commonly  used  windows,  the  variation  in  wkl  +  Wrkl-—) 
between  the  minimum  and  maximum  values  is  less  than  a  factor  of  two,  considerably 
less  than  the  several  orders  of  magnitude  variation  in  the  value  of  wkl  alone.  The  result 
is  a  substantial  decrease  in  the  noise  enhancement  effect. 


Some  examples  of  the  effect  on  SNR  for  different  windows  are  listed  in  Table  1.  In  this 
case  the  change  in  SNR  was  calculated  by  appropriately  modifying  (15)  to  get 


A  SNR 


—10  log 


(M 

2D  1 

M  t^(.W(Dk-D+l)+^(Dk-D+l+f))2 


(21) 
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Table  1:  Adjacent  channel  filter  response. 


Window 

Change  in  SNR  (dB) 

Triangular 

0 

Hamming 

0 

Hanning 

0 

Kaiser  / 3  =  5 

-0.1 

Kaiser  (3  =  10 

-2.7 

Kaiser  (3  =  15 

-6.2 

Gaussian  a  =  2.5 

-0.6 

Blackman 

-1.8 

Chebyshev  R  =  —70  dB 

-1.8 

B  ohman 

-2.2 

Nuttall 

-3.8 

Blackman-Harris 

-4.0 

Of  the  commonly  used  windows,  the  triangular  window  satisfies  (24)  exactly  and  the 
Hanning  and  Hamming  windows  nearly  satisfy  it  (the  variation  between  the  minimum 
and  maximum  values  of  wkl  +  uy/i:|  _  m  .  is  less  than  0.2%).  Adjusting  the  weights  for 
the  other  windows  to  satisfy  (24)  can  be  accomplished  using 


w'kl 


wkl 


Wki  +W{ki_M)  -  1 
2 


(22) 


,  ^+W{ki- M)-l 

W{kl-M)  =  ^(fcl-f) - 2 -  (23) 

which  adjusts  both  weights  equally.  This  adjustment  modifies  the  properties  of  the 
window,  however,  and  tends  to  make  the  window  properties  more  “hanning”-like. 


In  the  event  that  the  weights  satisfy  the  constraint 


=  1 


(24) 


the  noise  enhancement  effect  disappears  completely  and  the  estimation  of  xk  simplifies 
to 


Xk2 


zmn  +  {  l)'1  1) 

Dp™-1 


(25) 


The  filtering  approach  introduced  here  is  called  the  spectral  filter  throughout  the  rest  of 
this  report.  A  few  additional  comments  about  using  this  approach  arc  also  in  order.  As 
the  spectral  filter  has  been  presented  up  to  this  point,  the  decimation  rate  is  tied  to  the 
desired  filter  bandwidth,  Ffm„  according  to 


D 


M 

i‘2  —  h  +  1 


A 

Fbw 


for  i-2  >  i\.  The  output  data  sampling  rate  is  then  given  by 


F0  = 


Fs 

D 


(26) 


(27) 
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which  is  the  Nyquist  rate  (i.e.  F0  =  F}rw).  To  increase  the  output  sampling  rate,  while 
maintaining  the  same  filter  bandwidth,  rows  of  zeros  (zero  padding)  can  be  added  to  Z 
at  the  beginning,  end,  or  both.  The  resultant  new  decimation  rate  will  be  given  by 


D 


M  Fs 

- : - : -  <  — 

Tiz  +  *2  i  ]  T  1  FjjW 


(28) 


where  nz  is  the  number  of  rows  of  zeros  added  and  it  is  assumed  that  nz,  i\,  and  i 2  are 
chosen  to  yield  an  integer  value  for  D.  The  new  output  sampling  rate  is  still  calculated 
using  (27). 


2.3  Processing  Requirements 

In  this  section,  the  computational  requirements  of  the  proposed  spectral  filter  are 
considered:  specifically,  the  arithmetic  requirements.  Additional  processing  due  to 
housekeeping,  memory  transfers,  etc.,  are  not  considered  as  this  becomes  very 
dependent  on  the  processor  performing  the  filtering,  the  type  of  memory  used,  and  the 
way  the  algorithm  is  implemented  in  software,  which  is  beyond  the  scope  of  this  report. 
For  the  purposes  of  comparing  the  proposed  approach  with  other  approaches  discussed 
later  in  this  report,  it  can  be  assumed  that  a  generic  processor  is  used  in  all  cases  (i.e. 
not  optimized  in  a  way  which  favours  one  approach  over  another),  critical  memory 
transfers  involve  cache  memory  only,  and  that  all  approaches  are  efficiently  coded  in 
software. 

In  order  to  consider  the  processing  requirements  of  the  proposed  filtering  method,  it  is 
useful  to  define  a  generic  operation  which  can  be  used  to  quantify  the  computations 
required  to  carry  out  a  given  process.  To  do  this,  consider  the  following  summation 
series 

Q 

E  U(iV(i  ^ 

3=1 

which  is  representative  of  most  of  the  mathematical  computations  performed  in  the 
filtering  process.  There  are  three  cases  to  consider:  both  uq  and  vq  are  real-valued;  one 
is  real-valued  and  the  other  is  complex  valued;  and  both  are  complex-valued. 

In  the  first  case,  where  both  are  real-valued,  0  real  multiplications  and  Q  —  1  additions 
are  required  to  carry  out  (29).  Given  this  is  the  simplest  case,  it  is  useful  to  define  an 
operation  as  a  single  real  multiplication  plus  a  single  real  addition.  Hence,  assuming 
that  Q  >  1,  which  will  normally  be  the  case,  then  Q  operations  are  required  to  carry 
out  (29). 

In  the  second  case,  where,  for  example  uq,  is  real-valued  and  vq  is  complex  valued, 
then  (29)  can  be  represented  as 

Q  Q  Q 

E  uq{^{vq}  +  3^{vq})  =  E  Uq^{vq}  +  3  Y  Uq%{Vq}  (30) 
q= i  q= i  9=i 
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which  requires  twice  the  number  of  operations  as  the  first  case,  or  2  Q  operations. 

In  the  third  case,  where  both  uq  and  vq  arc  complex  valued,  then  (29)  can  be 
represented  as 

Q 

+  3%{Vq}) 

9=1 

Q  Q 

=  9f{«,}9f{vg} 

9=1  9=1 

Q  Q 

+  3  XI  ^{Uq}^{vq}  +  j^2$i{uq}$t{vq} 

9=1  9=1 

(31) 

which  requires  four  times  the  number  of  operations  as  the  first  case,  or  4 Q  operations. 

Finally,  some  of  the  filter  processing  may  also  include  Q  real  multiplications  without 
any  additions.  In  this  case  the  number  of  operations  is  considered  to  be  )qQ  for  this 
report.  This  is  based  on  the  example  of  the  Intel  P4  processor  which  can  carry  out 
floating  point  additions  with  a  throughput  of  2  clock  cycles  and  a  latency  of  4  clock 
cycles,  and  floating  point  multiplication  with  the  same  throughput  but  a  latency  of  6 
clock  cycles.  Under  ideal  conditions,  for  a  large  number  of  pipelined  instructions,  the 
average  processing  time  per  addition  will  then  be  the  same  as  the  average  processing 
time  per  multiplication.  In  real  systems  this  ideal  will  likely  not  be  achieved, 
particularly  when  the  performance  of  alternate  hardware  implementations  and  different 
software  libraries  is  taken  into  account.  Hence  the  assumed  equivalence  between 
addition  and  multiplication  times  is  only  an  estimate,  so  that  the  ’’number  of 
computations”  results  that  follow  should  also  be  considered  estimates.  However,  in 
terms  of  assessing  the  relative  performance  of  different  filtering  approaches  in  terms  of 
the  number  of  computations,  the  results  remain  valid. 

Having  now  examined  and  quantified  an  “operation”,  the  processing  requirements  of 
the  proposed  spectral  filtering  method  can  be  determined.  The  processing  can  be 
broken  into  four  separate  steps  which  can  be  examined  separately  to  determine  the 
individual  requirements.  These  four  steps  arc: 

1.  the  application  of  the  window  weighting  coefficients  in  (2), 

2.  the  application  of  the  FFT  on  the  columns  of  Y  to  get  S  shown  in  (6), 

3.  the  inverse  FFT  on  the  columns  of  Sm  to  get  Z  shown  in  (8),  and 

4.  the  generation  of  the  filtered  output  values  from  Z  using  (14). 

In  the  first  step,  the  application  of  the  weighting  coefficients  is  a  straightforward 
process  involving  the  multiplication  of  each  of  M  x  K  complex  data  values  by  a  real 
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weighting  coefficient.  Given  that  a  multiplication  involving  a  complex  and  a  real  value 
is  equivalent  to  two  real  multiplications,  and  remembering  the  factor  used  for 
multiplication  without  addition,  then 


C  i 


2  MK 
2 


2  N  -  M 


(32) 


where  (5)  was  used  to  replace  K. 

For  the  second  step,  the  computational  requirements  of  an  M -point  FFT  arc  given  by 

M  log  M  (33) 

complex  multiplies  and  adds  [2],  which  is  equivalent  to 


4 M  log  AM 


(34) 


real  multiplies  and  adds.  Hence,  the  number  of  operations  for  converting  the  K 
columns  of  Y  to  S  using  the  FFT  is  given  by 

C2  =  AKM  log2  M  =  (8N  —  AM)  log2  M  (35) 

where  (5)  was  used  to  replace  K. 

The  third  step  involves  an  inverse  FFT  carried  out  on  a  subset  of  the  rows  of  S,  namely 
Sm.  The  processing  requirements  are  given  by 

M  M  8 N  -  AM 

C3  =  AK  —  log2  —  =  - — - (log2  M  —  log2  D)  (36) 

where  the  inverse  FFT  size  is  M/D  and  the  relationship  expressed  in  (5)  was  used  to 
replace  K  (which  represents  the  number  of  columns  of  Sm). 


For  the  fourth  step,  the  number  of  filtered  output  values  L  which  are  generated  is  given 
by  the  number  of  unique  values  of  x\~  in  Z,  or 


L  = 


N 

~D 


(37) 


Assuming  that  processed  data  from  the  previous  and  next  blocks  of  Z  are  available, 
then  (20)  can  be  used  to  process  all  the  data  (i.e.  the  data  in  the  first  half  of  the  first 
column  of  Z  is  combined  with  data  in  the  last  column  of  the  previous  version  of  Z  and 
the  data  in  the  last  half  of  the  last  column  of  Z  is  combined  with  data  in  the  first  column 
of  the  next  version  of  Z  using  an  appropriately  modified  version  of  (20)  without 
resorting  to  (14)).  Noting  that  (20)  involves  a  complex  addition  and  multiplication 
(where  the  inverse  of  the  denominator  term  can  be  computed  beforehand),  then  the 
processing  requirements  are 


C4  =  AL 


AN 

~D 


(38) 
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In  terms  of  the  overall  processing,  the  first  two  processing  steps  arc  only  carried  out 
once,  regardless  of  the  number  of  signals  to  be  demodulated.  Since  these  two  steps  are 
also  being  carried  out  for  detection  purposes  as  well,  then  arguably  these  steps  could  be 
considered  as  “free”  processing  and  not  included  as  paid  of  the  filter  processing 
calculations.  However,  since  current  realtime  systems,  which  generate  time/frequency 
spectral  data  for  detection  purposes,  would  likely  use  contiguous  (not  overlapping)  data 
blocks  then  a  50%  overlap  would  result  in  a  doubling  of  the  processing  1 .  Hence,  the 
added  processing  attributable  to  the  filtering  process  will  be  \{C\  +  C2). 

The  last  two  processing  steps  must  be  carried  out  for  each  signal  to  be  demodulated. 
Consequently,  the  added  processing  requirements  due  to  these  two  steps  will  be 
71(63  +  C4)  where  n  is  the  number  of  signals  to  be  demodulated. 

From  this,  the  total  number  of  computations  required  for  the  proposed  filtering  method 
is  given  by 

C  =  i(Ci  +  C2)  +  n(C3  +  CA) 

=  i  (2 N  -  M  +  {8N  -  4 M)  log2  M) 

( 8N  —  4M  ^  4  N\ 

+n  (  - — - (log2  M  -  log2  D)  +  —  J 

I  zj.77 

=  -(2 N  -  M)(  1  +  41og2  M)  +  —  ((2 N  -  M)( log2  M  -  log2  D)  +  N) 

(39) 

For  most  filtering  applications  where  signal  demodulation  is  the  objective,  large  blocks 
of  input  data  will  be  used  so  that  it  is  reasonable  to  assume  N  3>  M.  This  leads  to  a 
simplification  in  the  expression  for  the  total  number  of  computations  according  to 


C 


N  +  4 N  log2  M  + 


4niV 

D 


(2  log2  M 


2  log2  D  +  1) 


(40) 


The  cost  of  computing  the  Fourier  coefficients  ej27rA:,/A/  and  e-o2nk/M  for 
k  =  0, M  —  1  has  not  been  considered  since  they  only  need  to  be  computed  or 
loaded  from  a  table  once  and  are  used  many  times  (i.e.  since  it  is  assumed  that  N  3>  M 
then  K  1)  so  that  the  additional  processing  load  will  be  negligible. 


'Note  that  overlapping  leads  to  improved  detection  performance  depending  on  the  FFT  window  used  [1], 


DRDC  Ottawa  TM  2005-198 


13 


3.  Temporal  Filtering 


In  the  following,  fast  variants  of  a  single  stage  filter  and  a  multiple  stage  filter  arc 
presented  which  arc  appropriate  for  demodulation  purposes.  Their  processing 
requirements  arc  also  discussed  in  a  way  which  makes  comparisons  with  the  spectral 
filter  discussed  in  Section  2  straightforward.  As  discussed  previously  in  Section  2.3,  the 
processing  requirements  arc  restricted  to  arithmetic  operations  only. 


3.1  Single  Stage  Filter 


The  basic  filter  operation  can  be  represented  by 

M 

zi  =  E  hkVk+i- 1  (4!) 

fc=l 


where  z\, ...,  zn-m+ i  represents  the  filtered  output  data,  hi, ...,  Hm  represents  the  low 
pass  filter  coefficients  (taps)  which  are  chosen  according  to  the  desired  filter  response, 
and  yi,...,yN  represents  the  input  data  after  the  signal-of-interest  has  been 
downconverted  to  baseband  according  to 


Vk 


(42) 


The  downconversion  is  done  prior  to  filtering  so  that  only  a  single  set  of  filter 
coefficients  is  required,  regardless  of  the  center  frequency,  fc,  of  the  signal.  This 
simplifies  storage  and  processing  requirements  as  well  as  ensuring  that  the  filter 
response  remains  constant. 

In  terms  of  the  processing  requirements,  the  downconversion  features  N  complex 
multiplications,  which  is  equivalent  to  4Ar  real  multiplications  and  2 N  real  additions,  or 

Ci  =  3  N  (43) 


operations.  It  also  involves  the  generation  of  the  mixing  term  represented  by 
e-j2irfc(k-i)/Fa '  por  qlc  pU1p0ses  0f  speed  it  is  assumed  that  the  appropriate  values  of 
the  mixing  terms  are  tabulated  in  cache  memory  2.  The  number  of  operations  will  then 
be  given  by 

N 

C2  =  -  (45) 

'This  can  be  done  by  saving  the  sequence 

a(n)  =  e~j2nnAf/F‘  (44) 

for  n  =  0, ...,  2^—1.  The  quantity  A /  is  the  frequency  quantization  increment  and  is  chosen  small  enough  that 
replacing  fc  by  mAf  does  not  significantly  degrade  the  accuracy  of  the  filtering  process.  The  mixing  sequence  is 
then  given  by  s(0),  s(m),  s(2m), s(mN  —  m)  where  the  fact  that  s(Fs/Af  +  fc)  =  s(k)  is  also  used.  The  main 
computational  requirement  is  the  generation  of  the  table  index  terms  0,  m,  2m, ...  which  require  either  a  single  real 
multiplication  or  singe  real  addition  per  term 
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In  terms  of  producing  the  filtered  output,  for  many  applications  it  is  not  necessary  to 
generate  every  possible  output  value  of  Zk,  but  rather  the  decimated  sequence 
represented  by  z\,z\+d,  Z1+2 d,  ■■■■  Hence,  calculating  only  the  required  output  values 
leads  to  the  modified  filter  equation  given  by 

M 

zi+jD  =  ^  hkUk+jD  (46) 

k= 1 

where  the  index  j  =  0, L  —  1,  and  z\+jjj  is  the  decimated  output.  In  this  case,  L 
output  values  arc  produced  so  that  the  number  of  operations  required  to  generate  the 
decimated  sequence  is 

C3  =  2  LM  (47) 

where  hi,...,  Hm  are  assumed  to  be  real-valued. 


The  total  number  of  operations  is  simply  the  sum  of  (43),  (45)  and  (47)  times  the 
number  of  signals,  n,  to  be  demodulated,  or 

7n  N 

C  =  n(C1  +  C2  +  C3)  =  — - b  2nLM  (48) 

An  alternate  form  of  this  expression  can  be  derived  by  considering  that  the  number  of 
output  data  values  will  be  given  by 


L 


1  + 


N  —  M 
D 


(49) 


To  simplify  later  comparisons  between  different  filtering  approaches,  it  is  assumed  that 
N  is  chosen  so  that 


L 


1  + 


N  —  M 
D 


(50) 


The  total  number  of  operations  can  then  be  rewritten  as 

^  7nN  n  ,  r  2nM(N  —  M) 

C  = - h  2  nM  + - ^ 

2  D 


(51) 


Finally,  using  the  assumption  that  N  S>  M,  this  simplifies  to 


C  R 

7nN 

2nM  N 

i  _  + 

2 

D  | 

(52) 


3.2  Multistage  Filtering 

3.2.1  Halfband  Filter 

Faster  filter  implementations  can  be  achieved  by  using  a  multistage  filter.  In 
particular,  a  simple  halfband  filter  can  be  used  to  successively  halve  the 
bandwidth  of  the  data  at  each  stage  until  the  desired  bandwidth  is  achieved. 
The  filter  equation  for  each  stage  is  given  by  modifying  (46)  to  get 

m2 

zi+jDk{k)  =  E  hm  ym+jDk(k )  (53) 

m=  1 
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for  k  =  1, K  where  K  is  the  number  of  filter  stages,  A)/,:  =  2k  is  the 
cumulative  decimation  of  the  data  carried  out  by  the  kth  stage  (i.e.  using  a 
halfband  filters,  each  stage  decimates  the  data  by  a  factor  of  two  so  that  the 
cumulative  decimation  by  the  kth  stage  is  2k),  M2  is  the  number  of  filter 
coefficient  for  each  halfband  filter,  and  y,ik)  and  Zi(k)  represent  the  input  and 
output  data  of  the  stage  filter,  respectively  (hence  for  k  >  1,  then 
Vi(k)  =  Zi{k  -  1)). 

Before  considering  the  computational  requirements,  it  is  first  useful  to 
determine  the  number  of  decimated  output  values  available  at  the  output  of 
each  stage.  Letting  L 1,  L2 ■ ....  Lk  represent  these  quantities,  then  (50)  can  be 
rewritten  as 


Lk  —  1  —  2 2^k~l  (^4) 

for  k  =  1, ...,  K,  where  represents  the  number  of  filter  output  values  for 
the  kth  stage  (as  well  as  the  number  of  input  values  for  the  next  stage),  and 
Lq  =  N  and  Lk  =  L.  Beginning  with  L\,  and  successively  solving  then 

Li  =  l-^M2  +  ^iV  (55) 

L2  =  1  —  -M2  +  -L\ 

=  l  -  +  \n  (56) 


L3  —  1  —  -M2  +  -L2 


7  7  1 

= - M2  +  -N 

4  8  8 


(57) 


and  so  on.  For  the  kth  filter  stage,  the  general  solution  is  given  by 

Lk  =  1  —  2^2  + 

2k+i  _  2  2k  -  1  1 

=  2 k  2^  M<2  +  ¥n 


=  (2  -M2)  + 


N  +  M2-2 

2k 


(58) 


for  k  =  1, ....  I\. 
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The  number  of  computations  is  found  by  taking  the  requirements  for 
downconverting  the  signal  given  by  (43)  and  adding  this  to  the  filter 
requirements  found  by  applying  (47)  to  each  stage  and  summing  the  results. 
This  is  then  multiplied  by  the  number  of  signals,  n,  to  be  demodulated 
yielding 


C 


7nN 

2 

7  ?iN 
2 

7  nN 
2 

7nN 

2 


K 

+  n  'y  ^  2L);M2 
k=  1 

,,  ^  /  ,,  JV  +  M2-2\ 

+  2raM2  E(2- m2  4  ^  ) 

fc=i  '  ' 

+  2 nM2K (2  -  M2)  +  2nM2(N  +  M2  - 
+  2 nM2K (2  -  M2)  +  2nM2{N  +  M2  - 


K  1 
2)E^ 


(59) 


where 

D  =  2X  (60) 

is  the  decimation  rate  for  all  K  stages. 


It  is  relatively  easy  to  show  that  the  multistage  filter  implementation  as 
discussed  so  far  actually  leads  to  a  slower  filter  implementation  than  the  single 
stage  filter  discussed  previously.  The  main  advantage  of  using  a  multistage 
filter  implementation  comes  from  taking  advantage  of  the  fact  that  almost  half 
of  the  coefficients  of  a  halfband  lowpass  FIR  filter  arc  zero  if  M2  is  odd.  For 
example,  an  eleven-tap  filter  designed  using  a  Hamming  window  has  the 
coefficients 


[/ii, ...,  hn]T 


0.00506 

0 

-0.0419 

0 

0.288 

0.497 

0.288 

0 

-0.0419 

0 

0.00506 


(61) 


where  every  second  coefficient  is  zero  with  the  exception  of  the  central 
coefficient.  From  this  it  may  be  concluded  that  for  M2  =  3,  7, 11, 15, ....  the 
number  of  nonzero  coefficients  is  given  by 


M' 


M2  +  3 
2 


(62) 
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It  is  not  necessary  to  consider  the  cases  M2  =  5,  9, 13, 17, ...  since  this  results 
in  h\  =  0  and  h\j2  =  0  making  the  filter  equivalent  to  one  with  M2  —  2 
coefficients. 

The  zero-value  coefficients  mean  that  fewer  multiplications  arc  required  at 
each  stage3  so  (59)  is  modified  to  become 

C  = 


where  it  as  not  necessary  to  modify  the  expression  for  L\~  since  it  is  based  on 
the  total  number  of  coefficients,  not  the  number  of  nonzero  coefficients. 

Since  N  3>  M2  for  most  applications,  then  the  number  of  operations 
simplifies  to 


7nN 

2 

7nN 


K 


+  n  'y  ^  2LfcM2 


k= 1 


2  +  2nM^K(2  -  M2)  +  2 nM'2(N  +  Af2  -  2)  f  1  -  — 


7nN 

2 

7nN 


+  2 n  [  K{2  -  M2)  +  2 n(  ( N  +  M2  -  2) 


2  ; 


2  ) 


^  —  nK(M2  +  3) {M2  —  2)  +  n{M2  +  3){N  +  M2  —  2)  (  1  —  — 


7r?  AT  (  1 

C  ~  —  +  nlV(M2  +  3)  M  -- 


(64) 


For  comparative  puiposes,  the  number  of  coefficients  for  each  halfband  filter 
can  be  related  back  to  the  number  of  coefficients  for  an  equivalent  single  stage 
filter  according  to 


M2  =  ^1  +  1  (65) 

which  follows  from  (58)  since  N  =  M  when  L  =  1.  It  is  assumed,  in  this 
case,  that  M  and  D  arc  chosen  so  that  this  expression  yields  a  value  for  M2 
from  the  set  {3,  7, 11, 15, ...}.  Finally,  substituting  this  expression  back  into 
(64),  yields 


C  s 

15niV 

(M-  5\ 

4  +  nJS 

2 

V  0  J 

(66) 


which  allows  direct  comparisons  between  the  multistage  filter  and  the  other 
approaches  previously  discussed. 

3Provided  the  software  implementation  takes  advantage  of  the  zero  values  -  some  popular  software  libraries  do 
not 
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Note  that  further  reductions  in  the  number  of  computations  arc  possible  by 
modifying  the  number  of  filter  taps  on  a  stage  by  stage  basis,  where  a  lower 
number  of  taps  is  used  in  the  earlier  stages.  This  modified  halfband  filter 
approach,  however,  requires  more  careful  design  and  is  not  considered  in  the 
rest  of  this  report. 

3.2.2  Composite  Filter 

One  of  the  drawbacks  of  the  multistage  halfband  filter  for  demodulation  is  that 
the  filter  bandwidth  is  tied  to  the  decimation  rate  according  to 

Fbw  =  Fs/D  (67) 

Hence,  if  the  signal  is  filtered  to  the  correct  bandwidth,  the  output  sampling 
rate  of  the  filter  will  equal  the  Nyquist  rate.  Many  digital  demodulation 
techniques  require  the  signal  to  be  oversampled  by  a  factor  of  two  or  more,  so 
this  presents  a  problem. 

One  possible  solution  is  to  use  fewer  filter  stages  to  achieve  the  desired  output 
sampling  rate.  However,  this  results  in  a  wider  filter  bandwidth  and  a  greater 
susceptibility  to  noise  and  cochannel  interference.  A  more  complete  solution 
is  to  reduce  the  number  of  stages  to  achieve  the  correct  sampling  rate  and  then 
follow  this  by  a  one  stage  non-decimating  filter  designed  to  reduce  the 
bandwidth  to  the  desired  value.  The  computational  penalty  of  adding  a  single 
stage  filter  is  not  high  since  the  filter  is  working  on  decimated  data,  not  the  full 
input. 

When  implemented  this  way,  the  appropriate  halfband  filter  size  is  still 
determined  using  (65)  except  with  the  modification 

M  -  1 

M2  =  - - — +  1  (68) 

Umax  -L 

where  Dmax  represents  the  maximum  decimation  rate  (at  any  higher 
decimation  rate  the  signal  will  be  undersampled)  and  is  given  by 

Dmax  =  p-  (69) 

“bw 

The  number  of  computations,  C\  due  to  the  multistage  halfband  filter  alone  is 
given  by  (64),  which  is  repeated  and  modified  here  as 

Ci  « 


7~'y~  +  nN(M2  +  3)  fl  -  ^ 


7nN 


+  nN 


M  -  1 
Fmax  1 


+  4 


(70) 
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where  D  is  the  actual  decimation  rate  (D  =  2h  where  K  is  the  number  of 
stages). 

In  assessing  the  computational  penalty  of  the  final  stage  filter,  both  the 
number  of  input  data  values  and  filter  size  arc  required.  The  number  of  input 
data  values  feeding  the  final  stage  filter  will  equal  the  number  of  data  values 
output  by  the  multistage  halfband  filter  and  is  found  using  (58)  where  k  =  K. 
Taking  into  account  the  latest  relationship  for  M2  expressed  in  (68),  then  (58) 
can  be  rewritten  as 


Lr 


1  + 


N  -1 
D 


M  -  If  D-  1 

D  \  Diri.n.x  1 


(71) 


which  simplifies  to 


Lr  ~ 


N 

~D 


for  large  N .  The  filter  size,  Mi,  can  also  be  derived  from  (71)  since 
L[{  =  Mi  when  N  =  M.  Accordingly, 


Mi 


1  + 


M  -  1 
D 


Dmax  D 
Dmax  1 


(72) 


(73) 


Using  the  derivation  of  (47)  as  a  guide,  and  given  that  for  large  N  the  number 
of  input  data  values  L  j<  will  approximately  equal  the  number  of  output  values 
(since  it  is  a  non-decimating  filter),  then  the  number  of  computations  required 
for  the  final  stage  is  given  by 


U2  ~  2nL  kM\ 


2 nN  |  2nN(M  —  l) 

1 T  +  W2 


Or. 


Dr 


(74) 


The  total  number  of  computations  of  the  composite  filter  will  be  the  sum  of 
the  contributions  from  each  of  the  filters,  or 

C  «  Cl  +  C2  (75) 

Combining  (70)  and  (74),  and  then  simplifying  yields  the  final  result 


15  nN  2nN  nN(M  —  1){D2  —  3D  +  2Dmax) 
2  D  D'2(Dmax  —  1) 


(76) 
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4.  Comparative  Results 


In  the  following  two  sections,  the  frequency  response  and  computational  requirements 
of  the  single  stage,  multiple  stage  (halfband  and  composite),  and  the  spectral  filters  arc 
compared  and  analyzed.  These  comparisons  arc  carried  out  with  emphasis  on  the 
demodulation  of  signals  which  arc  narrowband  relative  to  the  input  data  bandwidth. 

4.1  Frequency  Response 

The  filter  frequency  response,  of  an  ideal  filter  would  pass  all  signal  energy  within  a 
desired  passband,  while  completely  rejecting  all  interfering  signals  and  noise  outside 
this  band.  Attempting  to  achieve  this  kind  of  performance  comes  at  the  cost  of  large 
filter  lengths  as  illustrated  in  Figure  3.  In  this  example,  the  frequency  response  of  a 
polyphase  filter4  is  shown  for  different  filter  length.  To  make  the  results  as  general  as 
possible,  the  filter  length  is  defined  as 

M  -  /»(£)  (77, 

where  (3  is  the  filter  length  factor. 


Figure  3:  Filter  response  as  a  function  of  filter  length. 

The  relevant  features  of  the  filter  response  arc  the  flatness  of  the  passband  (the 
frequency  region  between  -0.5  and  0.5),  the  roll-off  rate  in  the  transition  band  (i.e.  the 
rapid  drop  in  filter  gain  immediately  outside),  and  the  height  of  the  sidelobes  in  the 
stopband.  A  flat  response  across  the  passband  minimizes  distortion  to  the  desired  signal 

4The  filter  was  designed  based  on  the  procedure  outlined  in  [1]  and  using  a  Hanning  window. 
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while  a  rapid  roll-off  and  low  sidelobes  minimizes  cochannel  interference  and  noise. 
Better  response  can  be  achieved  by  increasing  the  filter  length  as  illustrated  in  Figure  3. 
Flowever,  increasing  filter  length  increases  the  processing  requirements  (i.e.  C 
increases  as  a  function  of  M  in  (40),(52),and  (66)),  hence  the  smallest  filter  length  is 
usually  chosen  where  acceptable  performance  is  still  achieved. 

What  is  meant  by  “acceptable  performance”  will  depend  on  the  application  as  well  as 
the  signal  environment.  Specifying  a  maximum  frequency  width  to  the  transition  bands 
while  maintaining  reasonable  passband  flatness  and  stopband  attenuation  is  the  criteria 
often  used.  Given  the  trend  observed  in  Figure  3,  an  approximately  equivalent  criteria 
requires  a  given  attenuation  at  a  frequency  offset  of  ±0.75  x  Ff)li:  relative  to  the  filter 
center  frequency.  The  critical  points  are  marked  by  the  vertical  dotted  lines  in  Figure  3. 
In  Figure  4,  the  filter  gain  at  ±0.75  is  plotted  for  various  values  of  (3.  From  this  figure, 
and  choosing  30  dB  criteria  as  a  reasonable  attenuation,  acceptable  performance  is 
achieved  when  (3  =  8. 


Figure  4:  Filter  gain  at  a  normalized  frequency  of  ±0.75  as  a  function  of  filter  length. 

The  choice  of  (3  =  8  is  not  meant  to  be  considered  as  “the”  choice  since  various 
engineering  considerations  come  into  play  in  its  determination.  For  example,  the  choice 
of  criteria  parameters  is  as  much  based  on  the  signal  environment  as  on  user 
preferences.  A  larger  transition  band  requirement,  for  example,  would  lead  to  a  smaller 
filter  size.  Additionally,  there  are  some  filter  designs  which  generate  lower  sidelobes 
(hence  better  stopband  attenuation)  which  might  lead  to  a  smaller  filter  size.  The  price, 
however,  is  poorer  passband  response  and  a  larger  transition  band.  The  filter  design 
chosen  here  was  considered  to  be  a  reasonable  compromise. 

Using  (3  =  8,  and  given  Fs/ Ff)W  =  256,  examples  of  the  frequency  response  of  a  single 
stage  filter,  a  multistage  halfband  filter  (a  composite  filter  was  not  required  for  this 
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example),  and  the  spectral  filter  arc  shown  in  Figure  5.  The  best  performance  is 
achieved  using  the  single  stage  filter.  This  is  not  surprising  since  it  has  the  greatest 
flexibility  in  design  given  that  all  M  coefficients  of  the  filter  may  be  chosen 
independently.  The  least  flexibility  is  found  using  the  multistage  halfband  filter  since 
only  M2  filter  coefficients  arc  available.  By  comparison,  the  spectral  filter  has  M/2 
weighting  coefficients  which  can  be  independently  modified  if  the  constraint  in  (24)  is 
followed.  Performance-wise,  the  composite  filter  has  a  similar  roll-off  rate  to  the  single 
stage  filter  but  higher  sidelobes.  The  spectral  filter  has  a  better  response  with  lower 
sidelobes  but  a  slower  roll-off  than  the  single  stage  filter. 


(a) 


(b) 


(c) 


Figure  5:  Frequency  response  of  (a)  the  single  stage  filter,  (b)  the  composite  filter,  and  the  (c)  time/frequency 

spectral  filter.  The  filter  size  factor  was  (5  =  8. 
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4.2  Computational  Comparisons 


Before  comparing  the  computational  advantages  and  disadvantages  of  the  filters 
discussed  so  far,  it  is  useful  to  define  the  relationship  between  the  decimation  rate  and 
the  signal  bandwidth,  which  is, 


D 


1 

a 


(78) 


where  a  is  the  oversampling  factor.  For  digital  demodulation  techniques,  the  complex 
sampling  rate  of  the  filtered  signal  will  obviously  have  to  be  at  least  equal  to  the  signal 
bandwidth  (a  =  1).  However,  better  demodulation  performance  can  often  be  achieved 
using  higher  sampling  rates  (typically  a  >  2). 


Using  (3  =  8  and  a  =  2  as  reasonable  starting  values.  Figures  6  and  7  plot  the  effect  of 
varying  these  values  on  the  number  of  computations  for  the  three  filter  types  considered. 
The  number  of  computations  has  been  shown  normalized  with  respect  to  the  number  of 
input  samples.  Additionally,  it  is  also  assumed  that  only  one  signal  requires 
demodulation.  The  stepped  nature  of  the  curve  for  the  multi-stage  halfband  filter  is  due 
to  the  fact  that  the  halfband  filter  sizes  arc  restricted  to  M2  =  5, 9, 13, 17, ....  Hence  M2 
was  selected  by  first  using  (68)  and  then  selecting  the  closest  legal  value. 


Examining  the  results,  the  multistage  filter  requires  only  about  25%  of  the  number  of 
computations  of  the  signal  stage  filter  in  both  Figures  6  and  7.  The  reasons  will  be 
discussed  shortly. 

Comparing  the  single  stage  filter  to  the  spectral  filter,  the  single  stage  filter  required  a 
smaller  number  of  computations  for  lower  values  of  / 3  and  a,  but  a  greater  number  for 
higher  values.  This  is  due  to  the  fact  that  increasing  either  the  filter  size  or 
oversampling  factor  results  in  a  linear  increase  in  the  number  of  computations  for  the 
single  and  multistage  filters,  but  significantly  less  for  the  spectral  filter.  In  fact,  the  main 
computational  load  for  the  spectral  filter  is  generating  the  time/frequency  spectral  data 
(the  first  term  in  (40))  so  that  increasing  the  filter  size  results  in  a  logarithmic  increase 
in  the  computations  (i.e.  log  M)  and  increasing  the  oversampling  factor  (or  decreasing 
the  decimation  rate)  has  only  a  minor  effect. 


Comparing  the  multistage  filter  to  the  spectral  filter,  the  multistage  filter  required  fewer 
computations  in  all  cases  shown.  As  with  the  single  stage  filter,  there  is  a  linear  increase 
in  the  number  of  computations  as  a  function  of  filter  size.  However  in  Figure  6,  the 
filter  size  was  not  increased  high  enough  for  the  spectral  filter  to  obtain  better  results. 


The  computation  equations  representing  the  filters  under  evaluation  (i.e.  equations 
((40),  (52),  and  (76))),  can  be  approximated  by  simpler  forms  which  arc  more  suitable 
for  evaluation  if  certain  conditions  arc  met.  Particularly  if  the  desired  filter  bandwidth 
is  a  small  percentage  of  the  input  sampling  rate,  and  the  oversampling  factor  is  not  large 
so  that 


D2  >  Dmax  >  1 


(79) 
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5  10  15  20  25  30  35 

Filter  Length  Factor  (P) 

Figure  6:  Computation  ratio  C/N  as  a  function  of  the  filter  length  for  an  oversampling  factor  a  =  2  and  a  single 

signal. 


Oversampling  Factor  (a) 


Figure  7:  Computation  ratio  C/N  as  a  function  of  the  oversampling  factor  for  a  filter  length  factor  /3  =  8  and  a 

single  signal. 
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Additionally,  if  the  number  of  signals  to  be  filtered  is  also  not  large  (1)  n),  then  the 

simplified  computation  rates  arc  given  by 


cv 

N 


4  log  2  M 


(80) 


Css  7 n  |  2 nM 

W  ~  Y +  d 


(81) 


Col 

N 


15  n  nM 
~ 2  ^  ~n 


(82) 


which  represent  the  spectral,  single  stage,  and  composite  filters,  respectively.  For 
wideband  systems  processing  narrowband  signals,  these  assumptions  arc  reasonable. 


From  these  expressions,  it  is  apparent  that  the  lower  computation  rate  of  the  composite 
filter  compared  to  the  single  stage  filter  is  primarily  due  to  the  advantage  of  zero 
coefficients  in  the  halfband  filter  (note  the  extra  factor  of  2  in  the  second  term  of  (81) 
compared  to  the  second  term  of  (80))  and  the  beneficial  effect  of  doing  the  decimation 
as  early  in  the  filtering  process  as  possible  (another  D/Dmax  improvement). 

It  is  also  clear  from  these  expressions  that  (as  expected)  increasing  the  number  of 
signals,  n,  results  in  a  linear  increase  in  the  number  of  computations  for  the  single  stage 
and  composite  filters,  but  no  significant  increase  in  the  number  of  computations  for  the 
spectral  filter.  Re-examining  Figures  6  and  7  in  this  context,  then  for  two  or  more 
signals  the  spectral  filter  will  generally  be  faster  than  the  single  stage  filter  (e.g.  for  two 
signals  the  ratio  C/N  in  either  figure  will  be  doubled  for  all  single  stage  results  making 
them  greater,  for  the  most  paid,  than  the  corresponding  FFT-based  results  which  remain 
unchanged). 

Comparing  the  spectral  filter  to  the  composite  filter, with  respect  to  the  number  of 
signals,  is  slightly  more  complicated.  Mathematically,  the  spectral  filter  is  faster  when 

Ctf  <  Cef 


15n  nM 
4  log2  M  <  — +  — - 

^  D max 


Given  (69)  and  (77),  then  the  condition  can  be  rewritten  as 

log 2  M  < 


(83) 


(84) 


Figure  8  illustrates  this  condition  graphically  over  a  range  of  different  values  for  M,  (3, 
and  n  (note  that  in  generating  these  curves,  no  attempt  was  made  to  ensure  that  M2  was 
a  legal  value).  Although  these  are  approximate  results,  they  suggest  that  for  reasonable 
filter  sizes  (M  <  65536)  and  passband  response  (/?  >  8),  the  spectral  filter  is  faster 
when  four  or  more  signals  are  processed  simultaneously. 
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Figure  8:  Parameter  curves  (f3  versus  n)  which  mark  the  boundaries  between  faster  processing  for  the  spectral  filter 
(above  the  curve)  and  faster  processing  for  the  composite  filter  (below  the  curve).  Different  curves  are  plotted  for 

different  filter  sizes  (M ). 


5.  Conclusions 


In  this  report  a  fast  filtering  approach  suitable  for  applications  involving  the 
simultaneous  demodulation  of  a  number  of  narrowband  signals  captured  in  wideband 
digital  data  was  introduced.  This  filtering  approach  is  tailored  for  spectrum  surveillance 
systems  which  first  generate  a  time/frequency  spectrum  from  the  wideband  data  for  the 
purposes  of  signal  detection  and  tracking  purposes.  Using  this  spectral  data,  the  data  is 
then  processed  to  produce  the  desired  filtered  and  decimated  signals  at  baseband 
frequencies.  This  is  different  from  the  normal  approach  where,  for  example,  an  FIR 
filter  is  used  to  process  the  input  wideband  data  to  produce  the  desired  processed  signal. 
In  this  report,  the  new  approach  is  called  the  spectral  filter. 

The  main  advantage  of  the  spectral  filter  is  that  the  bulk  of  the  processing  involves  the 
generation  of  the  time/frequency  spectrum  which  can  be  done  efficiently  using  a  fast 
Fourier  transform  (FFT).  Subsequent  processing  to  extract  and  filter  the  desired  signals 
requires  very  little  additional  processing,  hence  there  is  relatively  little  difference 
whether  one  or  a  number  of  signals  arc  extracted.  By  way  of  comparison,  for  normal 
time  domain  based  filters  the  processing  increases  as  a  linear  function  of  the  number  of 
signals. 

The  main  disadvantage  is  that  each  time/frequency  spectrum  needs  to  be  generated 
from  blocks  of  data  which  overlap  the  previous  data  block  by  50%.  Although  this  also 
has  advantages  for  signal  detection,  for  the  sake  of  speed,  current  systems  such  as 
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MiDAS  use  no  overlap  so  that  spectral  filtering  requires  that  the  processing  be  doubled. 
For  the  time  domain  filters,  there  is  no  such  doubling  since  the  spectral  data  is  ignored. 

Two  other  filters  were  also  introduced  for  comparative  purpose,  namely,  a  single  stage 
FIR  filter  and  a  multistage  FIR  filter.  The  single  stage  filter  has  the  greatest  flexibility  in 
terms  of  designing  the  filter  frequency  response.  The  multistage  filter  has  considerably 
less  flexible,  but  is  computationally  more  efficient. 

In  terms  of  filter  frequency  response,  the  spectral  filter  falls  in  between  the  single  stage 
and  multistage  filters.  It  is  more  flexible  and  has  a  better  response  than  the  multistage 
filter,  but  is  still  relatively  inflexible  when  compared  to  the  single  stage  filter. 

In  terms  of  computational  speed,  both  the  single  stage  and  multistage  filters  are  faster 
when  only  one  signal  is  filtered.  However  when  two  signals  arc  processed 
simultaneously,  the  spectral  filter  is  faster  than  the  single  stage  filter,  and  when  four 
signals  arc  processed  simultaneously,  the  spectral  filter  is  also  faster  than  the  multistage 
filter. 

The  assessment  of  comparative  computational  speeds  is  based  on  the  assumption  that 
implementing  the  spectral  filter  requires  a  doubling  of  the  spectral  processing  already 
being  carried  out.  If  no  spectral  processing  was  previously  carried  out,  then  the  relative 
performance  will  be  worse,  and  the  number  of  signals  for  which  the  spectral  filter  is 
faster  becomes  four  and  eight  (instead  of  two  and  four  stated  in  the  previous 
paragraph).  On  the  other  hand,  if  the  full  spectral  processing  is  already  being  carried 
out  (i.e.  the  input  data  blocks  overlap  by  50%),  then  the  spectral  filter  will  be  faster  for 
any  number  of  signals. 

The  assessment  of  comparative  computational  speeds  is  also  based  on  ignoring  various 
non-arithmetic  processor  activities  such  as  memory  transfers,  pipelining,  etc..  These 
activities  could  alter  the  relative  performances  of  the  various  filtering  approaches 
analyzed.  In  addition  some  filter  approaches  lend  themselves  better  to  efficient 
hardware  designs  which  could  significantly  alter  the  relative  performance  results  from 
those  predicted  by  a  simple  number  of  computations  assessment.  Actual  testing  of  the 
various  filtering  approaches  on  real  systems  would  be  required  to  investigate  this. 

The  main  conclusion  is  that  for  typical  applications  where  the  demodulation  of  four  or 
more  narrowband  signals  is  desirable,  the  spectral  filter  is  the  preferred  choice. 
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